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Abstract 

A  number  of  practically  relevant  problems  involving  the  impulsive  motion  or  the  rapid 
rotation  of  bodies  immersed  in  fluid  are  susceptible  to  vortex-like  iiistability  modes. 
Depending  upon  the  configuration  of  any  particular  problem  the  stability  properties  of 
any  high-wavenumber  vortices  can  take  on  one  of  two  distinct  forms.  One  of  these  is  akin 
to  the  structure  of  Gortler  vortices  in  boundary  layer  flows  whilst  the  other  is  similar  to 
the  situation  for  classical  Taylor  vortices. 

Both  the  Gortler  and  Taylor  problems  have  been  extensively  studied  when  crossflow 
effects  are  excluded  from  the  underlying  ba»se  flows.  Recently,  studies  have  been  made 
concerning  the  influence  of  crossflow  on  Gortler  modes  and  here  we  use  a  linearised 
stability  analysis  to  examine  crossflow  properties  for  the  Taylor  mode.  This  work  allows 
us  to  identify  the  most  unstable  vortex  as  the  crossflow  component  increases  and  we  show 
how,  like  the  Gortler  case,  only  a  very  small  crossflow  component  is  required  in  order  to 
completely  stabilise  the  flow.  Our  investigation  forms  the  basis  for  an  extension  to  the 
nonlinear  problem  and  is  of  potential  applicability  to  a  range  of  pertinent  flows. 


^  Research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  contract 
No.  NASl-19480  while  the  author  was  'n  r'*'''dcncc  at  the  Institute  for  Computer  Applications  in  Science  and 
Engineering  (ICASE),  NASA  Langley  Research  C^enter,  Hampton,  VA  2.3681 
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1.  Introduction 

The  problem  of  predicting  vortex  breakdown  is  of  great  importance  within  the  field 
of  trajisition  research.  In  many  practical  situations  it  is  desirable  to  influence  breakdown 
characteristics:  for  example  in  aerospace  design  the  requirement  is  that  transition  be 
delayed  as  long  as  possible  in  order  to  reduce  drag  whereas  in  jet  intakes  turbulence 
is  promoted  in  order  to  improve  fuel  mixing.  In  numerous  applications  where  flow 
curvature  is  present  the  stability  problem  is  clo.sely  related  to  structures  describing  Taylor 
and/or  Gortler  modes  at  large  Taylor /Gortler  numbers.  Undoubtedly  in  many  cases  the 
t  xperimental  flow  is  fully  three  diuicn.sioxial  ux  character  and  until  recently  little  theoretical 
work  had  been  conducted  on  such  23roblems. 

The  linear  stability  of  short  wavelength  Taylor-Gortler  vortices  was  examined  by  Hall 
[1]  for  both  fully  developed  and  boundary  layer  flows.  The  model  jiroblem  considered  was 
the  flow  contained  between  a  pair  of  concentric  cylinders  whose  separation  is  much  smaller 
than  the  radii  of  the  cylinders.  A  WKB  type  analysis  shows  that  in  the  Dean  problem  (in 
which  both  cylinders  are  at  rest  and  the  fluid  is  driven  by  a  constant  azimuthal  i>ressure 
gradient)  short  wavelength  vortices  are  constrained  to  lie  in  a  thin  layer  located  at  a 
known  position  between  the  cylinders.  In  contrast,  for  the  Taylor  jiroblem  (in  which  the 
outer  cylinder  is  at  rest  and  the  inner  cylinder  rotates  uniformly)  the  vortex  structure  is 
traj323ed  at  the  inner  cylinder  wall.  These  two  idealised  cases  capture  the  essence  of  a  large 
number  of  im2)ortant  practical  flows.  For  exanq^le,  the  Dean  structure  forms  the  basis  for 
describing  the  evolution  of  Gortler  vortices  in  S2>atially  varying  boundary  layers  (although 
extra  non-23arallel  effects  do  need  to  be  considered)  whereas  the  wall-bounded  vortices  are 
relevant  to  a  number  of  impulsively  started  S23in-U23  2iroblems,  [2].  The  linear  stability 
theories  for  the  Taylor  and  Dean  cases  have  been  extended  to  weakly  and  fully  nonlinear 
accounts,  [3-6]. 

The  first  work  concerned  with  the  effect  of  adding  a  crossflow  com23onent  to  the  basic 
flow  was  conducted  by  Hall,  [7],  He  showed  that  only  a  very  small  crossflow  velocity  is 
needed  to  radically  alter  the  linear  stability  23roperties  of  short  wavelength  Gortler  vortices. 
In  many  23ractical  situations  interest  cannot  be  legitimately  restricted  to  high  wavenumljer 
modes  alone  and,  on  a  linear  stability  ba.sis,  it  is  frequently  the  fastest  growing  (or  most 
nnstaVile)  vortex  which  is  likely  to  l)e  important.  That  this  mode  is  of  interest  is  clear  for 
often  there  will  be  a  whole  S23ectrum  of  vortex  wavenumbers  that  may  be  excited  in  any 
one  2Jarticidar  23i'oblem  and  the  fastest  growing  of  these  is  the  one  that  is  dominant.  The 
most  unstable  Gortler  mode  within  a  two-dimensional  boundary  layer  v^as  id('iitifled  by 
Denier  et.  al.  [8]  and  the  linear  and  nonlinear  2)n)23erties  of  this  mode  as  the  crossflow 
com2'onent  is  incieasc’d  have  been  d«'velo2)ed  in  [9-11].  To  the  best  of  our  knowledge  th<’ 
most  unstable  vortex  in  the  Taylor  23robl<>m  has  yet  to  b('  not('fl. 

The  2)ur2)o.se  of  this  concise  article  is  to  outline  the  effect  of  crossflow  on  vortic('s  in 
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the  Taylor  situation  for  cylinders  with  a  small  .separation.  In  particular  we  shall  describe 
the  influence  of  crossflow  on  high  wavenumber  modes,  specify  the  structure  of  the  most 
unstable,  vortex  and  track  the  linear  stability  properties  of  this  mode  as  the  cros.sflow  grows. 
It  may  be  argued  that  this  problem  is  very  specific  and  therefore  is  of  little,  or  no,  practical 
use.  However  we  should  emphasise  that  the  analysis  here  encomijasses  the  ideas  behind 
a  large  number  of  important  cases  including  problems  involving  rapidly  rotating  bodies 
and  impulsive  motion.  Thus  the  merit  of  using  the  classical  small  gap  Taylor  problem  as 
a  paradigm  is  that  the  analysis  is  very  straightforward,  enables  the  important  features  of 
the  flow  development  to  be  seen  and  our  workings  may  be  easily  generalised  to  account  for 
more  complicated  situations. 

The  remainder  of  the  paper  is  organised  as  follows.  In  §2  we  describe  the  relevant 
stability  equations  and  in  §3  show  how  increasing  crossflow  modifies  the  large  wavenumber 
structures.  The  most  unstable  mode  is  obtained  in  §4  and  we  then  outline  the  influence  of 
crossflow  on  this  vortex.  We  close  with  a  brief  discussion. 

2.  The  governing  stability  equations 

The  linearised  disturbance  equations  pertaining  to  Taylor  vortices  have  been  derived 
many  times  and  so  here  we  restrict  ourselves  to  the  briefest  of  descriptions  of  these 
equations.  Consider  the  flow  of  incompressible  fluid  of  density  p  and  kinematic  viscosity  v 
in  the  gap  between  a  pair  of  concentric  cylinders  of  radii  R\  2ind  R2  {>  Hj).  We  assume 
that  the  curvature  of  the  gap  is  small;  that  is  6  =  (R2  —  Ri)/Ri  <C  1.  The  inner  cylinder 
rotates  with  angular  velocity  UofRi  whilst  the  outer  cylinder  is  at  rest.  It  is  convenient 
to  define  variables  x,y^z,t  by 

Ree'  r'-Ri  _  z'  Uot' 

""  S'  ^  R2-R1'  ^  R2-Ra'  Re{R2-R^y 

where  {r\d' ,z')  are  the  usual  cylindriceil  polar  co-ordinates,  t'  is  the  time  variable  and 
Re  =  Uq{R2  —  R\) I V  1  is  the  Reynolds  number.  Writing  the  velocity  field  (scaled 
on  the  characteristic  fluid  speed  Uq)  with  respect  to  (i,y,z)  as  {u,v,w)  and  scaling  the 
pressure  p  on  pU^  gives  that  in  the  basic  state 

(u,u)  =  (u,u)  =  (1  -  y,0)  +  O(^).  (2a) 

In  order  to  introduce  a  crossflow  component  into  the  flow  we  suppose  that  a  constant 
axial  pressure  gradient  acts  and  the  size  of  this  gradient  is  chosen  so  that  the  crossflow 
affects  the  subsequent  stability  equations.  It  is  a  routine  calculation  to  show  that  this 
occurs  when  the  axial  velocity  component  w  =  0(Re~^)  and  then  the  remainder  of  the 
basic  state  solution  is 

(w,p)  =  (w,p)  =  A  ^  {y(l  -  y)  +  0(S)}  {2z  -h  0(S)}"j  ;  (2b) 
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where  we  have  introduced  A  as  a  measure  of  the  strength  of  the  scaled  axial  flow.  The 
situation  discussed  in  this  article  has  been  studied  by  many  proliflc  authors,  the  reader  is 
referred  to  [12]  for  further  reading.  The  earliest  theoretical  and  experimental  studies  were 
apparently  made  by  [13]  and  [14]  respectively. 

The  Taylor  number  T  =  2Re^S  is  talcen  to  be  ari  0(1)  pareuneter  and  we  perturb  the 
basic  state  by  writing 


(u,  V,  w,p)  =  (u,  0,  iD,p)+A  |f7(r,  y,  z,  t),  ^Vix,  y,  2,  t),  ^W{x,  y,  2,  t),  V, 

(3) 


where  A  is  a  vanishingly  small  parameter.  If  we  substitute  (3)  into  the  linearised  forms  of 
the  continuity  emd  Navier-Stokes  equations  and  neglect  terms  of  relative  size  0{6),  we  ar  e 
left  with  the  four  disturbance  equations 


dU  dV  dW 

dx  ^  dy  ^  dz 


(4a) 


Vay2  dz^ 


U  =  a-yf£-V  +  Xy(l-y)§, 


\dy^  dz^ 


dP  dV  dV 


52  d'^ 


(1-y) 


dW 

dx 


+  A(l-2y)F  +  Ay(l-y) 


dW 
dz ' 


(46) 

(4c) 

(4d) 


which  need  to  be  solved  subject  to  the  no-slip  conditions  U  =  F  =  W  =  0ony  =  0, 1. 
In  the  absence  of  crossflow  (A  =  0)  arid  with  no  azimuthal  dependence  {d/dx  =  0)  system 
(4)  reduces  to  an  ordinary  differential  eigenvedue  problem  for  the  neutral  Taylor  number 
as  a  function  of  vortex  wavenumber  a.  It  is  well  known  that  for  a  <C  1  then  T  =  0(a~'^) 
whilst,  for  large  a,  T  =  0(0“*).  We  now  analyse  the  effect  of  crossflow  on  this  right-hand 
branch  of  the  neutral  stability  curve;  i.e.  on  neutral  large- wavenumber  vortex  structures. 


3.  The  effect  of  crossflow  on  large  wavenumber  modes 

It  was  shown  in  [1]  that  large  wavenumber  vortices  in  a  flow  with  zero  axial  conipt)nent 
are  trapped  in  a  thin  zone  of  depth  0(n“?)  next  to  the  inner  cylinder  wall  y  =  0. 
Analysis  of  eqiiations  (4)  reveals  that  when  the  scaled  crossflow  size  A  =  0(a)  this  wall 
structure  needs  to  be  adjusted  and  consequently  we  shall  examine  this  case  now.  In  the 
neighljourhood  of  the  neutral  stability  curve  it  is  known  that  the  vortices  have  temporal 
growth  rate  O(a^)  and  a  waviness  in  the  azimuthal  direction  on  an  0(a~3)  lengthscale. 
Using  the  scalings  proposed  in  [1]  we  write  the  Taylor  number 


r  =  roa"  +  r,a^  +  ..., 


(5a) 
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seek  modes  proportional  to 

El  =  exp  iaz  +  a^{(To  +  a~^ ai  +  . .  .)x  +  J  (fio(<)  +  -f  . .  ^  (5ft) 

and  allow  the  crossflow  size  A  =  aA.  If  we  further  define  the  0(1)  co-ordinate  ^  within  the 
wall  layer  by 

y  =  (6a) 

then  the  disturbance  quantities  {U,V,W,P)  develop  according  to 
{U,V,W,P)^ 

(jJo  +  a  kFo  o.Wi  . . . ,  Po  -|-  a^Pi  Ei . 

(6b) 

Substitution  of  (5)  &  (6)  into  (46,  c)  yields  at  leading  orders 


—(1  -f-  Q,o)Uo  =  —14, 


—(1  -|-  fio)14  —  ~PoUo, 


which  are  consistent  only  if 


To  =(l-f  fto)^ 


In  order  to  determine  the  structure  of  the  leading  order  disturbance  functions  it  is  necessary 
to  proceed  to  higher  orders  in  (4).  These  give  the  governing  equation 

(f2o  -t-  3)  —  Vo  ^  ^1  -f  fio  -l-  2iA^  —  ^^(1  +  ^o)  —  2(fi]  -I-  cto)^  =  0, 

which  is  a  scaled  version  of  the  Airy  equation  and  has  solution 

T/  A-  +  2xA'\  t  /'O  ^ 

(-  JTsh,  ') 


where 


(fio  +  3)  5 


2(fi]  -t-  (To) 


(l+fto  +  2zA)5  VTo . J 

hi  order  to  satisfy  the  inviscid  constraint  that  Vo  — 0  as  ^  >  0  it  is  necessary  to  choose 

C  so  that  Ai(—C)  =  0  whence 

C  w  2.3381  (8c) 

and  then  (86)  may  Vie  rearranged  to  give  the  Taylor  number  correction  Tj.  Although 
Vo  — ♦  0  as  ^  >  0  the  continuity  equation  (4a)  gives  Wo  =  i{dVo/dQ  and  so  a  passive 

viscous  layer  is  required  to  ensure  that  the  axial  disturbance  velocity  component  vanishes 
on  the  inner  cylinder.  The  details  of  this  layer  are  very  elementary  so  will  not  lie  pursued 
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Definition  (56)  shows  that  for  neutrally  stable  large  wavenumber  inodes  the  frequency 
terms  DojD]  are  necessarily  purely  imaginary.  Thence  result  (7)  gives  To  =  1  and  (86) 
shows  that 

Ti  =35C(l  +  4A2)3cos(^itan"^(2A))  .  (9) 

Consequently  for  large  scaled  crossflow  A,  T\  =  0(A3)  and  expansion  (5n)  implies  that 
when  A  ~  (z  the  neutral  Taylor  number  is  altered  by  an  O(a^)  amount  from  its  value  in 
the  zero  crossflow  state.  Furthermore,  the  extent  of  the  vortex  activity  is  compressed  into 
a  thinner,  0(a“’ ),  depth  region  adjacent  to  y  =  0  and  so  we  write 


A  =  a^A, 


n  ’x, 

(lOn, 6) 

^  J  (Do(f)  +  . .  ■)dt 

(10c) 

W  o  +  ...,a^.Po4-...^  E2 . 

(lOd) 

and  write  the  vortex  in  the  form 

iU,V,W,P)^  + 

Substitution  into  (4)  leads  to  the  disturbance  equations 


_  1  _  -  ^0  -  ilx)  -  l)  ^0  =  To^o, 

=  -Vo, 


which  need  to  be  solved  subject  toUo  =  Vo  —  dVo/dx  =  0  at  x  =  0  and  as  x  oo- 

At  this  point  it  is  worthwhile  to  recall  the  results  of  [7]  which  was  concerned  with 
describing  the  effect  of  crossflow  on  the  neutral  stability  properties  of  high  wavenumber 
Gortler  vortices.  In  the  absence  of  crossflow  in  that  case  the  mode  is  confined  to  a  thin 
layer  away  from  the  bounding  surface.  Three  distinctive  sizes  of  crossflow  were  identified. 
In  the  first  the  decay  of  the  vortices  away  from  the  thin  layer  is  primarily  due  to  the  fact 
that  the  local  Gortler  number  is  maximised  at  the  centre  of  the  layer;  at  a  higher  crossflow 
the  decay  is  driven  by  convective  effects  and  in  the  third  regime  the  layer  is  driven  to 
the  wall  and  the  vortex  is  governed  by  a  pair  of  viscous  equations  similar  in  style  to  (11) 
but  containing  quadratic  terms  as  well.  Thus  in  the  Gortler  (or  equivalently  the  Dean) 
problem  the  evolution  of  the  vortex  with  increasing  crossflow  passes  through  three  phases 
whereas  for  the  Taylor  problem  studied  here  we  have  only  two  stages. 

Equations  (11)  were  solved  using  ^  tinitp  difference  technique  and,  in  order  to  capture 
neutral  modes,  the  quantities  Do  and  do  were  taken  as  purely  imaginary.  The  results 
expressing  To  and  (Do  +  do)  as  functions  of  the  scaled  crossflow  A  are  shown  in  figure  1. 
Of  interest  is  the  solution  of  (11)  in  the  limits  of  large  or  small  cros.sflows.  As  A  — >  0  then 
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we  find  that  To  — >  1  and  Im(Oo  +  do)  — »  0  in  agreement  with  the  large  A  results  for  tlie 
ease  presented  above.  When  A  — >  ^  the  solution  ^livides  into  a  donbh  lioiindary  layer 
strnctnre.  We  postidate  that 

To  ~  d\  Hq  +  do  ~  icX  ( 12 ) 

for  some  real  constants  c  and  d.  In  the  majority  of  the  region  where  \  =  0(1)  then 

Uo  —  X  Lq  +  . . . ,  Vo  =  +  . . . ,  ( 13e ) 


where 


uci  =  v-„' 


and  the  solution  of  this  system  which  decays  exponentially  as  \  >  oc  is 

{ulv,h={x-'^KA\)^i\-^KAx))-  (14) 

where  is  the  modified  Bessel  fvmction  of  order  u  -  \Jd  +  \  (see  [15]).  As  \  — +  0  these 

^  —  i 

solutions  develop  singidarities  so  that  in  a  wall  layer  where  X  =  ^  =  (^(l)-  write 


Uo  =  \  C/+(r^  ,  'U-  \  ^  . . 


-  i  + 


(15e,/') 


These  functions  satisfy  the  coupled  equations 


\dY^  J  dY^  ”  \dY^  ) 


( 15c, d) 


+  +  + 

which  need  to  be  solved  subject  to  the  conditions  that  Uq  =  =  dV^  /dY  =  0  on  1'  =  0 

together  with  suitable  matching  conditions  to  (14)  a.s  F  — >  oo.  Numerical  solution  of  this 
problem  reveals  that 

c  = -4.706,  d  =  0.721,  (16) 


and  the  asymptotic  forms  (12)  for  To  and  Qq  +  do  are  sketched  on  figure  1.  As  may  be 
seen,  agreement  between  the  numerical  .solution  of  (11)  and  these  asymptotes  is  eminently 
satisfactory  and  demonstrates  that  as  the  (scaled)  crossflow  moves  through  this  O(a^) 
regime  the  changes  induced  in  the  neutral  Taylor  number  are  substantial.  Therefore  only 
a  very  small  axial  field  of  size  0{Re~^a^)  is  required  in  order  to  radically  distort  the 
neutral  curve  for  large  wavenumber  modes  and  oo  suggests  that  any  such  wavenumljer 
vortex  would  be  completely  stabilised  by  a  tiny  crossflow. 

In  many  practical  situations  Taylor/  Gortler  numbers  are  large  and,  whilst  the  effect  of 
any  crossflow  component  on  the  neutral  curve  is  of  interest,  it  is  also  important  to  describe 
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how  the  crossflow  iiiodifles  the  growth  rate  of  the  most  unstable  vortex.  As  was  mentioned 
in  the  introduction,  the  most  unstable  mode  is  frequently  the  one  that  is  dominant  in 
practice  and,  as  far  as  we  are  aware,  the  location  of  this  mode  has  not  been  writt<'n  down 
for  the  Taylor  problem.  Therefore  in  the  coming  section  we  first  obtain  the  most  unstable 
vortex  in  the  alisence  of  crossflow  and  then  consider  how  the  imposition  of  a  gradually 
increasing  crossflow  component  modifies  the  stability  properties  of  this  disturbance. 

4.  The  most  unstable  mode 

In  order  to  deduce  the  most  unstable  linear  vortex  in  the  al)sence  of  crossflow  it  is 
easiest  to  proceed  in  the  manner  of  Denier  ct.  al.  [8].  In  this  work  the  authors  considert'd 
both  the  growth  rate  of  inviscid  (0(1)  wavenumber)  Gdrtler  modes  as  the  wavenumber 
increases  and  also  the  properties  of  mode.s  in  the  vicinity  of  the  right  hand  branch  of  the 
linear  neutral  curve.  It  was  then  possilde  to  identify  an  intermediate  wavenumber  regime 
in  which  infinitesimal  high  Gdrtler  number  vortices  are  the  fastest  growing. 

The  inviscid  problem  relevant  to  Taylor  vortices  for  T  ^  1  has  been  considc'ri'd 
numerous  times  (see  [16]  for  example)  and  for  vortex  wavenumber  a  =  0(1)  tlu- 
distiirbance  takes  tb^  form 
(T\  F,  W,P)  = 

({/o(y)  +  ...,  r" Fo(y)  +  ....  THko(y)  +  ... ,  TPoiy)  +  . . .) exp  jnn-  +  7'^  J 

Sid).stitution  of  these  forms  into  the  governing  equations  (4)  shows  that  Fj  satisfies 

=  -«"(1  -  y)Vo.  (17) 

subject  to  the  boundary  conditions  that  Vq  =  0  at  y  =  0, 1.  This  gives  an  eigenproblem 
for  the  growth  rate  as  a  function  of  a.  In  order  to  locate  the  most  unstable  vortex  across 
the  entire  wavenumber  spectrum  it  is  useful  to  consider  the  l)ehaviour  of  the  solution  to 
(17)  as  a  — >  oo.  Elementary  analysis  (see  [2])  shows  that  the  eigenfunction  Iq  F  confined 
to  a  thin  zone  of  depth  0{a~^)  attached  to  y  =  0.  Thence 

D  =  1  -  |Ca-3 (18) 

where  the  constant  C  sati.sfles  Ai{  —  C)  =  0  and  has  been  given  in  (8r).  Thus  as  a  — >  oo 
the  growth  rate  of  the  vortex  asymptotes 

T'^  (l-\Ca-l  +  .  (19) 

Otto  [2]  ol)tainefl  this  result  for  his  problem  concerning  the  impulsive  s[)in-up  of  a  cylinder 
within  an  expanse  of  fluid  and  we  conclude  that  for  large  a  the  growth  rate  be<-omes 
independent  of  wavenumber  to  leading  order. 


8 


S.  R.  Otto  and  Andrew  P.  Bassorn 


R^'sults  (5a)  and  (7)  deinonstrate  that  at  large  vortex  wavenumbers  the  right  hand 
luaneh  of  the  neutral  stability  rurve  is  given  by  T  =  o'*  +  ....  Then  for  T  1  if  we 
examine  a  vortex  of  wavenumber  n  =  hT*  +  ...  with  a  =  0(1)  it  is  straightforward  to 
deduce  that  the  vortex  is  confined  to  a  region  of  depth  0{T~^  )  next  to  the  inner  cylinder 
wall.  Furthermore,  the  sizes  of  the  dist\irbance  components  {U,\\W,  P)  are  in  the  ratios 
(l,T5,r^,T5)  and  the  leading  order  vortex  growth  rate  term  is  0{T'^).  If  we  let  this 
growth  rate  0,  =  T^Cl  +  .  .  .  then  the  compatibility  criterion  arising  from  the  leading  order 
Ijalances  from  momentum  equations  (46,  e)  reduces  to 

h=\~cP. 

Thus,  as  woidd  be  expected,  the  vortices  are  unstable  for  d  <  1  and  stable  for  d  >  1.  In 
addition,  the  overall  growth  rate  is 

T5(l-d2)f...  (20) 

and  thus  the  expressions  (19)  and  (20)  become  comparable  when 

a  ~  i  . 

This  suggests  that  the  most  unstable  Taylor  vortex  has  wavenumber  )  and 

we  shall  confirm  this  presently.  Meanwhile  we  note  how  the  geometry  of  growth  rate/ 
wavenumber  sj)ac<*  for  the  Taylor  luoblem  contrasts  with  that  for  the  Giirtler  (or  Dean) 
case.  We  have  seen  from  (19),  (20)  that  for  a  large  wavenumber  range,  specifically 
1  ‘C  a  -C  T7,  the  growth  rate  of  the  Taylor  mode  is  constant  at  leading  order  and  only 
varies  in  higher  order  terms.  Conversely,  in  the  Gdrtler  case  the  leading  order  growth  rate 
increases  both  as  a  — >  cx:  in  the  inviscid  problem  and  as  d  — »  0  in  the  right  hand  branch 
cas‘-.  Therefore  in  the  interme(hate  regime  where  the  growth  rate  is  greatest  the  maximum 
is  a  much  more  pronounced  feature.  This  ditference  in  growth  rate  ))ehaviours  for  the  two 
prol)lems  does  have  important  consequences  for  th<‘  subsequent  analysis.  For  the  Gortler 
case  the  calculations  of  Bassoni  cC  Hall  [9]  show  that  steadily  incn’asiiig  crossflow  stabilises 
the  vortic<'s  and  only  a  small  crossflow  is  required  in  order  to  ensure  that  the  flow  is  stable 
to  infinitesimal  modes  irresja-efive  of  their  wavelength.  In  addition  it  was  found  that 
as  the  crossflow  grows  so  th<’  location  of  the  most  imstable  vortex  remains  in  the  same 
wavenumber  re’gime:  behaviour  which  is  not  repeated  for  the*  Taylor  cas<'  discus.sefl  below. 

In  order  to  prove  that  it  is  the  0(Tt%)  wavenumber  regime  which  coutains  the  most 
unstable  Taylor  mode  we  write 

(I  =  doT^,  (21a) 

and  seek  disturbances  of  the  form 

iu.  V.  w.P)  =  (to  + ....  tH''o  + ....  lio  +  ...,  r?  Po  +  .. .)  P,. 


(216) 
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wher 


;e 


£3  =  exp  iaoT^'^z +  T'^  J  (hoit)  +  T  5^1,0 +  •••) 


(It 


(21r 


<U(  i  t-'o  1  •  •  •  functions  of  7/  =  T»y.  Therefore  these  vortices  are  trapi)e(l  in  an  0{  T  >' 


>i.e  adjacent  to  y  =  0  and  we  shall  impose  a  crossflow  of  scaled  size 


A  =  TTo\. 


221 


This  scaling  is  arrived  at  Ijy  choosing  sufficient  crossflow  so  as  to  ijerturl;  the  vortex 
structure  relevant  to  the  zero  crossflow  problem  without  gros'-ly  changing  its  underlying 
characteristics.  We  substitute  (21),  (22)  into  equations  (4)  and  compare  like  powers  of  T. 
The  consistency  condition  arising  from  momentum  l>alances  (46,  c)  leads  to 


Clo  =  1, 


(23) 


as  would  be  anticipated  from  result  (20).  At  next  order  we  find  that  Vq  satisfies  the  scaled 
Airy  equation 


dr/ 2 


~  do(l  +  27 Ado) 


7/  + 


2(dQ  +  17] ) 
(It  27  Ado ) 


Demanding  that  Vo  =  0  at  r/  =  0  and  as  7/  ^  00  leads  to  an  expression  for  Q]  and  thence 
we  deduce  the  overall  vortex  growth  rate  Gr  =  T?  +  Cl\T^  +  o{T^ )  where 


C’(l  +  4d2A2)i 


2do 


|cos^|tan  '(2doA)^  +  7  sin  tan  ’(2doA)j  |  —  dg.  (24) 


In  the  absence  of  the  crossflow  component  Re{tt\ )  =  —  ^Cdp  ®  —  dg  which  clearly  — v  — oc 
as  uo  — >  0  or  as  do  00.  The  maximum  value  of  this  expression  is  — 4(C/6)5  ~  —1.973 
which  occurs  when  do  =  (C/6)5  0.702  and  consequently  we  have  the  result  that  in  the 

classical  Taylor  problem  the  the  most  unstable  mode  has  wavenumber  ~  0.702Tt^.  As  the 
scaled  crossflow  increases  then  Re{Cl\)  is  greatest  when  do  =  0(A“^)  and  |f2il  =  0(A3). 
Furthermore,  for  a  fixed  do,  as  A  increa.ses  so  the  overall  vortex  growth  rate  diminishes  aiifl 
crossflow  is  seen  to  have  a  stabilising  effect.  It  is  clear  that  a  revised  solution  structure 
is  needed  when  the  second  7crm  in  Gr  becomes  comparable  with  the  first,  i.e.  when 
17]  =  0(T5)  or  A  =  0(7^). 

We  therefore  define 

A  =  tU,  (25n) 

and  seek  the  wavenumber  of  the  most  unstable  vortex  for  this  case.  We  find  that  for  a  large 
range  of  wavenumbers,  specifically  1  n  T<,  the  leading  order  component  of  the  normal 
velocity  field,  say  Vo,  is  given  by  the  .solution  of  a  Whittaker  type  equation  (see  [15]).  Tin’ 
disturbance  structure  is  confined  to  lie  within  a  thin  zone  adjacent  to  y  =  0  whose  jirecise 
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depth  is  a  fnnction  of  th "  ’.vavenumber.  However  for  all  vortex  wavenumbers  1  <C  n  -C 
the  leading  order  disturbance  forms  are  determined  by  the  solution  of  a  common  equation 
which  may  be  solved  analytically  in  terms  of  the  modified  Bessel  function  A'^.  In  order  to 
satisfy  both  the  no-slip  conditions  at  the  wall  and  suitable  decay  conditions  at  the  edge  of 
the  thin  zone  it  is  found  that  the  leading  order  growth  rate  of  the  vortices  is  wnere 

fio  a  sohition  of  the  equation 

K^,  =0,  ^  i  +  A  ^  (256) 


The  dependeiice  of  flo  upon  A  is  shown  in  Figure  2.  When  A  =  0  vortex  modes  with 
wavenumber  1  <C  e  have  growth  rate  Ho  —  1-  As  A  increases  so  the  growth 

rate  diminishes.  Indeed,  as  A  I/n/2  so  Ho  — >  0  —  i/\/2  and  all  the  modes  with 
wavenumber  within  1  <C  u  -C  T?  are  made  neutral.  Again  therefore,  just  as  in  the  Gortl^r 
problem,  crossflow  is  seen  to  have  a  stabilising  influence  and  at  crossflows  of  greater  than 
/\/2  the  flow  configuration  is  stable  to  vortices  of  any  wavenumber  <C  T?. 

One  question  W(  have  not  addressed  here  is  that  concerning  the  structure  of  these  vortex 
modes  a-  Hie  crossflow  j^arameter  A  increases  through  the  critical  value  l/\/2.  Indeed 
when  A  —  l/\/2  precisely  then  the  formal  solution  lor  V'o  develops  a  singularity  and  the 
disturbance  structure  requires  a  critical  layer  type  zone  in  order  to  .smooth  this  out.  This 
situation  has  similarities  with  the  .scenaric  discussed  by  Blackaby  Choudhari  [17]  who 
were  concerned  with  inviscid  vortex  modes  in  slightly  three-dimensional  boundary  layers 
and  how  such  p<’rturbations  are  related  to  ♦’he  Taylor-Goldstein  equation  for  stratified  shear 
flows.  These  authors  found  solutions  for  which  the  leading  order  'vertical'  (i.e  normal  to 
wall)  velocity  component  of  the  inviscid  disturbance  becomes  infinite  at  the  critical  layer 
and  they  asserted  that  theirs  was  the  first  ca.se  in  which  this  type  of  behaviour  has  been 
found.  We  also  believe  this  to  be  the  case  and  suggest  that  we  have  here  another  problem 
where  this  phenomenon  is  seen. 

The  discussion  above  implies  that  for  a  large  range  of  vortex  wavenumbers  all  modes 
have  the  same  leading  order  growth  rate.  However,  such  vortices  are  not  the  fastest  growing 
ones  as  further  investigation  reveals  that  within  the  0(7’^)  wavenumber  regime  we  ha^'e 
the  unique  most  unstable  mode.  In  order  to  demonstrate  this  we  define  the  crossflow 
parameter  A  as  in  (25a)  and  look  for  disturbances  of  wavenumber  a  =  aT^ .  If  we  let 


E4  =  exp 


iaT^z  +  T-^ 


5  J'  (H(f)  +  ...)dt|, 


(26a) 


and  seek  solutions  of  the  type 

{U,V,W,P)r.{U +  ...,T^Va-...,T^W  +  ...,T^P  +  ..  .)E,, 


(266) 
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where  C  ,  W,  P  are  functions  of  the  sealed  co-ordinate  <p  =  T  ■<  y  we  find  that  U  and  V 
are  given  by  the  coupled  viscous  equations 


subject  to  the  conditions  that  U  —  V  ~  dV /d(p  =  0  at  =  0  aiifl  as  ci  — >  oo.  As  woidd  Ije 
expected  in  view  of  the  scalings,  these  equations  are  essentially  identical  to  those  relevant 
to  determining  the  right-hand  branch  of  the  neutral  curve,  (11). 

We  solved  (27)  for  the  (in  general  complex-  valued)  growth  term  Q  as  a  function  of  the 
crossflow  parameter  A  and  scaled  wavenumber  n  and  the  results  of  these  calculations  are 
summarised  in  figure  3.  This  figure  shows  that  for  A  <  1/V2  the  vortex  is  unstable  over 
a  finite  range  of  wavenuniljers  including  those  with  d  — v  0  and  the  extent  of  this  unstal.)le 
range  decreases  with  increasing  A.  Indeed  when  A  reaches  the  value  l/\/2  the  mode  with 
d  =  J  is  made  neutral  in  agreement  with  the  results  of  the  analysis  relevant  to  modes 
with  wavenumber  1  <C  «  <C  For  increasing  A  beyond  l/\/2  the  unstabh^  ))and  reduces 
further  until  when  A  %  0.942  the  whole  weivenuml)er  spectrum  is  made  stable. 

The  resvdts  given  in  figure  3  also  show  the  dependence  of  the  most  unstal)le  vortex 
wav«>nuniber  upon  the  crossflow  A.  If  we  denote  the  wavenumber  of  this  most  unstable 
mode  by  d„  then  as  A  — >  0  so  d„  0  since  we  know  fi'om  our  previous  workings  that  as 
the  crossflow  diminishes  so  the  most  unstable  vortex  reverts  to  within  the  0(TTc  )  regime. 
In  addition,  n,,  is  not  a  monotonic  function  of  A  and  as  A  — 0.942  it  is  the  mode  with 
wavenumljer  a  %  0.103  which  is  the  last  to  be  made  stable. 

The  calculations  in  this  and  the  preceding  section  give  related  but  not  identical 
information.  The  work  of  §3  may  be  regarded  as  describing  the  stability  })roperties  of 
a  vortex  of  specified,  large  wavenumber.  In  particular,  the  solution  of  equations  (11) 
accounts  for  the  dependence  of  the  neutral  Taylor  number  on  the  size  of  th»'  imposed 
crossflow  whereas  the  work  of  this  s«'ction  descriljes  how  crossflow  affects  the  stability  of 
vortices  within  a  flow  at  a  specified  large  Taylor  number.  We  notice  how  the  most  unstal)l(' 

3 

vortex  has  wavenumber  OlTfe  )  in  the  absence  of  crossflow;  with  increasing  crossflow  this 
fastest  growing  mode  moves  into  the  0{T^)  wavenumber  regime.  At  this  stage  th<'  right 
hand  neutral  branch  of  the  stability  curve  is  also  within  the  O(T^)  range  so  that  w«'  can  sts' 
that  as  the  crossflow  comixment  increases  from  zero  .so  tlie  location  of  the  most  unstal)l<' 
linearised  mode  moves  out  into  the  asymptotic  regim«>  containing  the  right  hand  I)ranch. 
This  contrasts  with  the  situation  in  the  Gdrth'r  ca.se  where,  for  a  large  Gdrtler  number  G 
and  in  the  absence  of  crossflow  the  right  hand  l)ranch  is  at  an  0(G5  )  wavenumber  whilst 
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the  fastest  growing  mode  has  wavenumber  0(^5),  [8].  As  the  crossflow  increases  the 
most  unstable  vortex  remains  within  the  0(G5 )  regime  and  the  right  hand  neutral  branch 
moves  into  this  regime  so  that  it  is  a  mode  of  wavenumber  0{G^ )  which  is  the  last  to  be 
stabilised. 


5.  Conclusions 

In  this  work  we  have  studied  how  the  imposition  of  a  small  crossflow  velocity 
component  aflFects  the  linear  stability  characteristics  of  vortices  in  high  Taylor  number 
flow.  In  particular  we  have  identified  the  most  unstable  vortex,  illustrated  that  for 
large  wavenumbers  a  crossflow  of  .size  0{Re~^  o^)  is  sufficient  to  completely  stabilise 
the  mode  and  proved  that  crossflow  of  makes  the  flow  immune  to  unstable 

vortices  irrespective  of  their  wavenumber.  In  some  ways  it  is  this  last  result  which  is 
the  more  significant  for  it  shows  that  in  many  flows  where  an  appreciable  degree  of 
three-dimensionality  is  j)resent  the  Taylor  vortex  mechanism  is  likely  to  l^e  unimportant 
compared  with  other  instability  forms. 

Our  results  may  i)e  compared  with  those  of  the  equivalent  Gortler  or  Dean  jjroblem 
in  the  following  way.  Suppose  that  instead  of  having  the  inner  cylinder  of  our  geometry 
moving  with  speed  Uq  it  is  at  rest  and  an  azimuthal  ]>ressure  gradient  is  imposed  of 
sufficient  strength  so  as  to  induce  an  0(Uo)  fluid  .speed  across  0  <  y  <  1 .  The  analysis 
of  [8]  suggests  that  in  the  absence  of  cros.sflow  the  most  unsta])le  vortex  has  wavenumber 
0(Ts ).  In  addition  a  crossflow  of  size  0{Re~^T^  )  would  Ije  needed  to  completely  stabilise 
the  flow-  a  value  which  is  seen  to  lie  greater  than  that  in  the  Taylor  problem  considered 
here.  This  suggests  that  in  situations  which  are  composed  partly  of  a  Taylor  type  situation 
and  partly  a  Dean  type  then  it  the  latter  component  which  is  dominant  in  determining  the 
staliility  characteristics  when  crossflow  is  pre.sent.  However,  there  will  still  lie  circumstances 
in  which  the  structure  presented  in  this  paper  is  important  and  such  cases  include  spin-uj) 
configurations  of  the  type  studied  liy  Otto  [2].  In  that  prolilem  a  cylinder  at  rest  in  an 
incompressilile  fluid  is  given  an  imimlsive  angidar  velocity.  A  Rayleigh  layer  is  set  uji 
and  this  flow  is  prone  to  vortices  which  are  identical  in  character  to  those  examined  here. 
If  a  liasic  flow  component  is  added  along  the  axis  of  the  cylinder  (so  this  could  model 
a  spinning  cylindrical  object  moving  through  a  fluul)  this  axial  flow  may  be  determim'd 
using  the  methofl  of  Glauert  Lighthill.  [18]. 

This  apiilicability  of  our  analysis  to  other  flow  jiroblems  is  an  attraction  of  our  working. 
The  flow  configuration  chosen  at  tlie  out.set  was  deliberately  taken  to  l)e  very  simple  l)ut 
it  does  cai)ture  most  of  the  essential  features  of  a  wide  class  of  flows  over  rotating  l)odies 
whenever  the  inviscid  Rayleigh  criterion  for  instability  is  most  violated  at  the  Ixxly  s\uface. 
Not  only  can  our  analysis  be  extended  to  other  geometries  Imt  it  may  also  lie  easily  arlapted 
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to  allow  for  vortex  motions  modulated  in  the  azimuthal  direction.  Consequently  it  has  the 
potential  of  proving  to  be  useful  in  a  number  of  related  problems. 

There  are  several  directions  in  which  our  calculations  might  be  extended.  The  most 
obvious  of  these  is  the  question  of  nonlinearity.  In  order  to  give  a  definitive  account  of 
crossflow  effects  it  would  be  necessary  to  allow  for  other  than  infinitesimal  vortex  modes 
by  pursuing  weakly  nonlinear  or  stiongly  nonlinear  computations.  Equivalent  work  for 
the  Gbrtler  problem  is  reported  in  [10,  11,  19]  where  it  is  shown  that  weak  nonlinearity 
staliilises  near  neutral  but  linearly  growing  modes  to  give  non-zero  finite  amplitude  vortex 
states.  Inclusion  of  full  nonlinearity  tends  to  suggest  that  the  vortices  break  down  in  a 
finite  distance  singularity  which  would  be  the  analogue  of  a  finite  time  singiilarity  in  the 
present  work.  In  addition,  of  theoretical  interest  would  be  a  rigf)rous  analysis  of  the  critical 
layer  structure  which  arises  just  as  vortices  of  wavenumber  1  <C  a  near  neutrality  as 

the  crossflow  approaches  Re~^T^ /\/2-  This  analy.sis  would  certainly  be  required  in  order 
to  investigate  nonlinear  configurations  near  this  cut-off.  However  the  practical  significance 
of  such  considerations  is  far  from  clear  at  this  juncture  for  at  this  crossflow  size  we  have 
shown  that  there  is  still  a  wide  range  of  higher  wavenumber  vortices  which  are  linearly 
unstable.  Indeed  it  is  not  until  that  the  axial  flow  becomes  greater  than  0.942/2e~’T5 
that  infinitesimal  Taylor  vortices  of  arbitrary  wavenumber  are  made  stable. 

We  are  grateful  to  Dr.  J.P.  Denier  for  his  valuable  comments  on  an  earlier  draft  of  this 
work.  The  research  of  the  first  author  was  supported  by  the  National  Aeronautics  and 
Space  Administration  under  NASA  contract  No.  NASl-19480  while  he  was  in  residence 
at  the  Institute  for  Computer  Applications  in  Science  and  Engineering  (ICASE),  NASA 
Langley  Research  Center,  Hampton,  VA  23681,  USA. 

References 

1.  P.  Hall,  J.  Fluid  Meek.  124  (1982)  475-494. 

2.  S.R.  Otto,  IMA  J.  Appl.  Maths.,  in  press  (1993). 

3.  P.  Hall,  J.  Inst.  Maths.  Applies.  29  (1982)  173-196. 

4.  P.  Hall  and  W.D.  Lakin,  Proc.  R.  Soc.  Land.  A  415  (1988)  421-444. 

5.  A.P.  Bassom  and  P.  Hall,  Stud.  appl.  Math.  81  (1989)  185-219. 

6.  J.P.  Denier,  IMA  J.  Appl.  Maths.  49  (1992)  15-33. 

7.  P.  Hall,  Proc.  R.  Soc.  Lond.  A  399  (1985)  135-152. 

8.  J.P.  Denier,  P.  Hall  and  S.O.  Seddougui,  Phil.  Trans.  R.  Soc.  Lond.  A  335  (1991) 
51-85. 

9.  A.P.  Bassom  and  P.  Hall,  J.  Fluid  Mech.  232  (1991)  647-680. 

10.  A.P.  Bassom  and  S.R.  Otto,  J.  Fluid  Mech.  249  (1993)  597-618. 

11.  S.R.  Otto  and  A.P.  Bassom,  Stud.  appl.  Math.,  in  press  (1994). 


14 


S.  R.  Otto  and  Andrew  P.  Bassom 


12.  R.  C.  DiPrima  and  A.  Pridor,  Proc.  R.  Soc.  Land.  A  366  (1979)  555-573. 

13.  S.  Goldstein,  Ptoc.  Camb.  Phil.  Soc.  33  (1937)  41-61. 

14.  R.  .1.  Corni.sh,  Proc.  R.  Soc.  Lond.  A  140  (1933)  227-240. 

15.  M.  Abramowitz  and  LA.  Stegun,  Handbook  of  Mathematical  Functions.,  (Dover,  New 
York  1972). 

16.  P.G.  Drazin  and  W.H.  Reid,  Hydrodynamic  Stability,  (CUP,  Cambridge  1981). 

17.  N.D.  Blackaby  and  M.  Choudhari,  Proc.  R.  Soc.  Lond.  A  440  (1993)  701-710. 

18.  M.B.  Glauert  and  M..J.  Lighthill,  Proc.  R.  Soc.  Lond.  A  230  (1955)  188-203. 

19.  .J.P.  Denier  and  P.  Hall,  J.  Fluid  Mech.  247  (1993)  1-16. 


‘  location  of  the  right-  hand 
’)  Im(fio  +  ^o)-  The  dotted 


Figure  3.  Sohitiou  of  viscous  equations  (27)  to  determine  leading  order  vortex  growth 
rate  Re{Cl)T^  as  a  fimction  of  vortex  wavenumber  aT*  and  crossflow  XT^ .  The  graphs 
correspond  to  crassflow  values  A  =  0.02,  0.05,  0.1  and  thereafter  at  intervals  of  0.1  up  to 
and  including  A  =  1.1.  The  A  symbols  indicate  the  locus  of  the  most  unstable  vortex 
waveniunber. 
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